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We study the nature of gauge fixing ambiguities in two dimensional gauge theories. We find that these ambi- 
guities can be related to the asociated spin model. They can be eliminated by means of a (multigrid) annealing 
algorithm. 



Gauge fixing is by now a useful tool in lat- 
tice gauge calculations. Smooth gauges, like the 
Coulomb gauge and the Landau gauge, are com- 
monly used in propagator calculations [ [p , as well 
as in wavefunction calculations [ ^] . In monopole 
studies, the so-called maximally Abelian gauge is 
frequently employed [|3j. However, since the work 
of Gribov [ H , it is known that 'fixing the gauge 
is not quite fixing the gauge'. In other words, it 
is possible to locally fix the gauge, but it need not 
always be possible to globally fix the gauge. Since 
one uses a slightly stronger condition on the lat- 
tice than in the continuum, it was long hoped that 
the Gribov ambiguity was absent on the lattice. 
A few years ago, it was shown that these hopes 
were idle [ |) . In order to get a better understand- 
ing of the nature of the Gribov ambiguity on the 
lattice, we shall study the occurrence of copies in 
2-d U(l) and 2-d SU(2) when fixing to the Landau 
gauge. 

The continuum version of the Landau gauge 
condition, 

0^»=O, Al = g\d, + AJg, 

translates to the lattice as 

Im J2 K(x) - U° f (x - A)]traceless = », (1) 

U9(x)=g^x)U tl (x)g(x + fL). 

As was mentioned above, usually a slightly 
stronger condition is used, namely minimizing a 
functional F, where F is given by 



(2) 



for SU(N) in d dimensions on a lattice containing 
V lattice sites. It is easy to see that if (||) is 
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Figure 1. The search for Gribov copies. R.G.T. 
stands for random gauge transformations 

minimized for a given gauge configuration U, (0) 
is automatically satisfied. 

The Gribov problem, finding multiple solutions 
for the Landau gauge fixing condition 

+ [D^A) , {dtf)g*\ = 0, 

Dfj,(A) being the covariant derivative, is equiva- 
lent to finding local minima for the functional F. 

The search for Gribov copies on the lattice is 
done as follows. We start from an initial con- 
figuration C\, and fix that configuration to the 
Landau gauge, keeping track of the gauge trans- 
formations involved. The resulting configuration 
is called L±. We then perform a random gauge 
transformation on C\, to obtain C2, and fix C2 
to the Landau gauge L2. We can repeat this as 
many times as we want. Since we have kept track 
of the gauge transformations gi between C\ and 
Li, we can construct the gauge transformation gij 
between Li and Lj-. 

9n{x)=g\{x)g j {x). (3) 
This procedure is displayed in figure [l]. 
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One can look upon the gauge fixing procedure 
from another point of view. Minimizing the func- 
tional F can be seen as finding the ground state of 
a spin system with pseudo-random couplings U, 
depending on both position and direction, with 
an action given by 



\ \ N N ~ 



Su(g) = -l3 s F(g,U). 



(4) 



We have put 'pseudo-random couplings U\ since 
only for gauge field coupling constant /3g = are 
the couplings U truly random. For finite /3g, the 
U are constrained to form plaquettes Up, which 
are to be taken from an equilibrium ensemble. 

Making a random gauge transformation is, in 
some sense, creating a spin configuration at spin 
coupling constant /3 S = 0. The standard pro- 
cedure to find the minimum of (Q) is employing 
a relaxation algorithm. Albeit impractical, this 
can be 'simulated' by performing Monte Carlo at 
fts — oo. We should therefore not be surprised if 
defects get trapped in the final, gauge fixed con- 
figuration: we don't expect a perfect crystal to be 
formed by dropping liquid iron into liquid helium. 

To determine whether we have reached conver- 
gence, we would like to use F. However, for non- 
trivial gauge configurations we do not know the 
the minimal value(s) of F. Therefore, we define 

fi — Fi-x — Fi. 

Fi is the value of F after i iterations. We can 
assign a relaxation time r to the minimization 
process by defining r by 



fi = Cexp(-i/r). 



(5) 



For 2-d U(l), the asociated spin model is the 
well known XY-model. It possesses a Kosterlitz- 
Thouless phase transition at a finite coupling (3 C . 
The system is described by vortices and antivor- 
tices for couplings (3 < (3 C , and by spin waves for 
couplings (3 > (3 C . 

We created pure gauge configurations by tak- 
ing Ufj,(x) = g'{x)g{x + ft), and fixed these to 
the Landau gauge. We used a standard relax- 
ation algorithm, combined with multigrid [|[. We 
experienced, even on a 16 2 lattice, difficulties in 
reaching the true Landau gauge, U^(x) = 1. On 
larger lattices it became virtually impossible to 
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Figure 2. The final, Landau gauge fixed configura- 
tion, starting from pure gauge on a 32 2 lattice. A few, 
widely seperated vortices and antivortices, as well as 
two spin waves, are left. 

reach the true Landau gauge. Figure || displays a 
typical gauge fixed configuration for a pure gauge 
configuration on a 32 2 lattice. 

Although the inclusion of overrelaxation in the 
minimization algorithm improves on the relax- 
ation time r and is able to remove all the vortex- 
antivortex pairs on the smaller lattices, the com- 
bined algorithm is not able to eradicate all pairs 
on the larger lattices. Nevertheless, in the follow- 
ing we shall make use of the combined algorithm. 

Turning to gauge fields taken out of an equilib- 
rium ensemble, we find basically the same results. 
We found that the r values depended slightly 
on the /3 value of the gauge fields, as well as 
on the lattice size, but not on the number of 
(anti)vortices or spin waves left in the gauge fixed 
configuration. 

For 2-d SU(2), the asociated spin model is 
the 0(4) non-linear sigma model. This model 
is asymptotically free and does not contain 
any topological objects such as instantons or 
(anti)vortices. Therefore, we do not expect to find 
any obstructions in fixing pure gauge to the Lan- 
dau gauge. This is confirmed by figure [| 

Taking gauge field configurations from an equi- 
librium ensemble, and fixing these to the Landau 
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Figure 3. 10 log F versus the number of iterations 
for pure gauge SU(2) configurations. An overrelaxed 
minimization algorithm was used in combination with 
multigrid with W-cycles. 
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Figure 4. 10 log / versus the number of iterations for 
one SU(2) gauge configuration. 

gauge, we at times do find Gribov copies. 

These copies behave differently from the ones 
we found for 2-d U(l). The different copies re- 
lax with different autocorrelation times t to their 
minima. This is depicted in figure |J. 

At fixed lattice size, it is observed that the re- 
laxation time t on average becomes larger when 
Pg is increased. Furthermore, the relative occur- 
rence of these copies becomes smaller; see table 

I 

We would like to find an observable / for 
the gauge transformation g^ to disinguish be- 
tween the different Gribov copies. This observ- 
able should be invariant under global spin rota- 
tions, i. e. f(gij) = f(gijh) , h € SU(2), since 
global spin rotations correspond to global gauge 
transformations U^x) — > U^ix) — h 1( Uf l (x)h. 
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Table 1. The number of times a gauge configuration 
yielded # copies for several combinations of the lattice 
size N and the gauge coupling f3 g . The /3 g values are 
chosen to keep the ratio N/^ a roughly constant. £ CT is 
the string tension correlation length. 

The first observable to satisfy this criterion is 
the most obvious one: the (spin model) energy 
density 

e(x) = Tr W( x )9( x + A) + .9 T (> ~ A)flOs)]- 

We use the parametrization 
g(x) = [.94l -iff ■ g] (x), 

since it allows us to change conveniently from the 
0(4) or vector notation, to the SU(2) or matrix 
notation, and vise versa. The second observable 
is given by 

Jfr) = ]Tj; 2 0r), 

i 

ji{x) = [e f _ l „e tj kigjd f _ l g k d L ,gi} (x). 

We have plotted 1 — \e{x) and J{x) for a gauge 
transformation taking one copy to the absolute 
minimum; see figure [| Surprisingly, the two fig- 
ures are almost identical in shape! The reason for 
this similarity is not (yet) understood and is being 
studied. 

The reason for copies to occur is that once we 
have come close to a local minimum, there is no 
chance to get away from this minimum, due to 
the minimization algorithm. This immediately 
implies a way out: using an annealing algorithm. 
The implementation for gauge fixing is clear from 
: we perform Monte Carlo on the spins g using 
the action Su{g) at some value of S . We should 
start with a (very) low value of 0s, and gradually 
increase (3s by small amounts A(3 S - After increas- 
ing (3 S , we should let the system reach equilibrium 
before increasing /3 S once more. When /3 S is large 
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Figure 5. J (top) and 1 — \e (bottom) for a 
gauge transformation relating a copy to the true Lan- 
dau gauge, obtained from a gauge configuration taken 
from an equilibruim ensemble at (3 = 48.0 on a 32 2 
lattice. 

enough, we can turn to ordinary relaxation, with 
or without overrelaxation. The starting value of 
Ps and A/3 S , as well as the value of [3 S where we 
can turn to the normal minimization must be de- 
termined in practice, and may depend on (3g. 

The interpretation of this process is straightfor- 
ward. Making a random gauge transformation is 
(in some sense) taking a spin configuration at in- 
finite temperature ((3 = T" 1 ). Performing Monte 
Carlo at low f3 s is nothing more than letting the 
(spin) system evolve in time, deep into the high 
temperature phase. Increasing /3 S is then lowering 
the temperature, which should be done gradually 
especially near a phase transition. If we have been 
careful enough, we can throw the system into liq- 
uid helium once we are deep into the low temper- 
ature phase. 

We used this annealing algorithm, using a 
multigrid algorithm as described in [ ^) for the 
Monte Carlo part, for both U(l) and SU(2) gauge 
fields. After tuning the aforementioned param- 
eters, we found no copies, even for very rough 
gauge fields; j3 = 24.0 on a 128 2 lattice for SU(2), 
where £ CT ~ 4. 

We can draw several conclusions from this find- 



ing. Firstly, the copies for SU(2) are probably a 
relic from the spin model. 

Secondly, the anology with spin glasses, which 
sometimes is made to explain gauge fixing ambi- 
guities, may be misleading. To be sure, it strictly 
holds for /3g = only, since for (3g ^ the cou- 
plings U are constrained, as was mentioned be- 
fore. Especially for larger values of (3g, we should 
expect that the local minima thin out, and be- 
come well seperated. This is indeed observed, but 
surprisingly, for fig values in the range of fig = 6.0 
- 192.0 on 32 2 lattices, we find AF^/F^ to be 
in the order of a few percent, i 7 ^ denotes the final 
value of F for a given configuration. This suggests 
that the deviation from a spin glass starts at lower 
values of /3g, and therefore, that the problem of 
finding the true Landau gauge is not as severe as 
finding the ground state of a spin glass. This fact 
encourages us to apply annealing algorithms to 
4-d gauge fixing problems. 
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